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Abstract — The Grassmannian space Q q {n,k) is the set of all 
k dimensional subspaces of the vector space F™. Recently, codes 
in the Grassmannian have found an application in network coding. 
The main goal of this paper is to present efficient enumerative 
encoding and decoding techniques for the Grassmannian. These 
coding techniques are based on two different orders for the Grass- 
mannian induced by different representations of /, '-dimensional 
subspaces of F™. One enumerative coding method is based on 
a Ferrers diagram representation and on an order for Q q (n, k) 
based on this representation. The complexity of this enumerative 
coding is 0(fc 5//2 (n — fc) 5//2 ) digit operations. Another order of the 
Grassmannian is based on a combination of an identifying vector 
and a reduced row echelon form representation of subspaces. The 
complexity of the enumerative coding, based on this order, is 
0(nk(n — k) log n log log n) digits operations. A combination of 
txfjthe two methods reduces the complexity on average by a constant 
^ factor. 

Index Terms — enumerative coding, Grassmannian, identifying 
qq vector, Ferrers diagram, partitions, reduced row echelon form. 

I. Introduction 

Let F q be a finite field of size q. The Grassmannian space 
(Grassmannian, in short), denoted by Q q (n, fc), is the set of 
all fc-dimensional subspaces of the vector space F™, for any 
'given two integers k and n, < fc < n. It is well known Q 
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that \Q q {n,k)\ = 
coefficient, defined by 
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is a g-ary Gaussian 
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= if k > n or k < 0. 
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, Coding (and related designs) in the Grassmannian was 
considered in the last forty years, e.g. [2-8|. Koetter and 
Kschischang |9j presented an application of error-correcting 
codes in Q q (n, k) to random network coding. This application 
has motivated extensive work in the area Q0-20|. A natu- 
ral question is how to encode/decode the subspaces in the 
Grassmannian in an efficient way. By encoding we mean a 
transformation of an information word into a fc-dimensional 
subspace. Decoding is the inverse transformation of the fc- 
dimensional subspace into the information word. 

To solve this coding problem, we will use the general enu- 
merative coding method which was presented by Cover ||2T|| . 
Let {0, l} n denote the set of all binary vectors of length n. 
Let 5 be a subset of {0, 1}™. Denote by ns{x\,X2, ■ ■ ■ ,Xk) 
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the number of elements of S for which the first k coordinates 
are given by (x±,X2, ■ ■ ■ ,Xk), where x\ is the most significant 
bit. A lexicographic order of S is defined as follows. We say 
that for x, y € {0, 1}™, x < y, if Xk < yk for the least index k 
such that Xk ^ yk- For example, 00101 < 00110. 

Theorem 1: OTTl The lexicographic index (decoding) of 
x G S is given by 



ind s (a;) = ■ n s (xi,x 2 , 



,0). 



Let 5 be a given subset and let i be a given index. The 
following algorithm finds the unique element x of the subset 
S such that inds(:E) = i (encoding). 

Inverse algorithm [21 1: For k — l,...,n, if i > 
ns(xi, X2, ■ ■ ■ , Xk-i, 0) then set Xk = 1 and i = i — 
ns(xi,x 2l . . . , ifc-i, 0); otherwise set Xk — 0. 

Remark 1: The coding algorithms of Cover are efficient if 
ns(xi,X2, ■ ■ ■ , Xj-i, 0) can be calculated efficiently. 

Cover [21] also presented the extension of these results to 
arbitrary finite alphabets. For our purpose this extension is more 
relevant as we will see in the sequel. The formula for calcu- 
lating the lexicographic index of x 6 S C {1, 2, 3, ... , M} n is 
given as follows. 



ind s (j;) = n s (xi. 



x 2 , 



,m). 



(2) 



j=l m<Xj 



Enumerative coding has various applications and it was 
considered in many papers, e.g. Ir22"1424l . Our goal in this 
paper is to apply this scheme to the set of all subspaces in 
a Grassmannian, using different lexicographic orders. These 
lexicographic orders are based on different representations of 
subspaces. Lexicographic orders also have other applications, 
e.g. in constructions of lexicographic codes (lexicodes) Il25ll . 

The rest of this paper is organized as follows. In Section ITfl 
we discuss different representations of subspaces in the Grass- 
mannian. We define the reduced row echelon form of a fc- 
dimensional subspace and its Ferrers diagram. These two con- 
cepts combined with the identifying vector of a subspace ifTSl 
will be our main tools for the representation of subspaces. 
We also define and discuss some type of partitions which 
have an important role in our exposition. In Section [III] we 
present a new lexicographic order for the Grassmannian based 
on a representation of a subspace by its identifying vector and 
its reduced row echelon form. For this order we describe an 
enumerative coding method, whose computation complexity is 
0(nk(n — k) log n log log n) digit operations per subspace. In 
Section [IV] we discuss the more intuitive order for the Grass- 
mannian based on Ferrers diagram representation and present 
a second enumerative coding method for the Grassmannian. 
In Section [V] we show how we can combine the two coding 
methods mentioned above to find a more efficient enumerative 
coding for the Grassmannian. In Section [VT] we summarize our 
results and discuss some related problems. 



2 



II. Representation of Subspaces and Partitions 

In this section we give the definitions for two concepts 
which are useful in describing a subspace in Q q (n,k): Ferrers 
diagram (which is defined in connection to a partition) and 
reduced row echelon form. Based on these concepts we present 
two representations for subspaces from which our enumerative 
coding techniques will be induced. Representation of subspaces 
is also important in other problems related to the Grassmannian. 
For example, in constructing error-correcting codes in the 
Grassmannian 111 812611 . 

A partition of a positive integer to is a representation of 
to as a sum of positive integers, not necessarily distinct. We 
order this collection of integers in a decreasing order. The 
partition function p(m) is the number of different partitions 
of to 1 11271281 . 

A Ferrers diagram T represents a partition as a pattern of 
dots with the i-th row having the same number of dots as the 
i-th term in the partition [ 1 27 28] (In the sequel, a dot will be 
denoted by a " • ")■ A Ferrers diagram satisfies the following 
conditions. 

• The number of dots in a row is at most the number of 
dots in the previous row. 

« All the dots are shifted to the right of the diagram. 

Remark 2: Our definition of Ferrers diagram (see |[T8l ) is 
slightly different form the usual definition 1 11271281 . where the 
dots in each row are shifted to the left of the diagram. 

A fc-dimensional subspace X g Q q {n, k) can be represented 
by a k x n matrix, whose rows form a basis for X. Such a 
k x n matrix is in reduced row echelon form (RREF in short) 
if the following conditions are satisfied. 

• The leading coefficient (pivot) of a row is always to the 
right of the leading coefficient of the previous row. 

• All leading coefficients are ones. 

• Every leading coefficient is the only nonzero entry in its 
column. 

For a given subspace X, there is exactly one matrix in RREF 
and it will be denoted by RE(X). For simplicity, we will 
assume that the entries in RE(X) are taken from Z q instead of 
F g , using an appropriate bijection. 

The Ferrers tableaux form of a subspace X, denoted by 
F(X), is obtained by removing from each row of RE(X) 
the leading coefficient and the zeroes to the left of it. All 
the remaining entries are shifted to the right. F{X) defines a 
unique representation of X. The Ferrers diagram of X, denoted 
by Tx, is obtained from F{X) by replacing the entries of 
F(X) with dots. 

Example 1: We consider a 3-dimensional subspace X QfV 7 2 
with the following 3x7 matrix in RREF given by 



1 1 1 
RE(X) = (0010101 
1 1 1 



Its Ferrers tableaux form and Ferrers diagram are given by 

110 • • • • 

J-{X) = 10 1 and Fx = • • • , respectively. 



Let \JF\ denote the size of a Ferrers diagram J 7 , i.e., the 
number of dots in T . A Ferrers diagram of a fc-dimensional 
subspace has size at most k ■ (n— k). It can be embedded in a 



k X (n — k) box. Let p(k, rj, to) be the number of partitions of 
to whose Ferrers diagram can be embedded into a box of size 
k x r). The following result was given in ll27l pp. 33-34]. 

Lemma 1: p(k, rj, to) satisfies the following recurrence rela- 
tion: 

p(k, r), to) = p(k, 7} — 1, to — k) + p(k — l,r), to) (3) 
with the initial conditions 
p[k, r), to) = if TO<0orTO>7y-/c and p{k, rj, 0) = 1. 

Let J 7 be a Ferrers diagram of size to embedded in a k x [n— 
k) box. The number of fc-dimensional subspaces whose Ferrers 
diagram is F, is equal to q m . By (Q]) this implies the following 
theorem fl] p. 327] which shows the connection between the 
g-ary Gaussian coefficients and partitions. 

Theorem 2: For any given integers k and n, < k < n, 
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k(n—k) 
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where a m — p(k, n — k, to). 

The order defined in Section [IV] is based on Theorem [2] We 
order the subspaces by the size of their Ferrers diagrams. The 
order of Ferrers diagrams with the same size is explained in 
Section [IV] Two subspaces with the same Ferrers diagrams are 
ordered lexicographically by their Ferrers tableaux forms. This 
order seems to be the most natural order of G q (n, k). But, a less 
natural representation, which follows, and its related order, will 
lead to a more efficient enumerative coding. 

Each fc-dimensional subspace X 6 Q q (n,k) has an identify- 
ing vector v(X) ifTHl . v(X) is a binary vector of length n and 
weight fc, where the ones in v(X) are exactly in the positions 
(columns) where RE(X) has the leading coefficients (of the 
rows). 

Let X £ G q {n, k) be a fc-dimensional subspace. The ex- 
tended representation, EXT(X), of X is a (fc + 1) x n 
matrix obtained by combining the identifying vector v(X) = 
(v(X) n , . . and the RREF RE{X) = (X n , . . .,Xx), 

as follows 



EXT(X) = 



v(X), 
X n 



v(X) 2 

x 2 



v{X) x 



column vector 



Note, that v(X) n is the most significant bit of v(X). Also, Xj 
is a column vector and v{X)i is the most significant bit of the 

( " { S- )■ 

Example 2: Consider the 3-dimensional subspace X of Ex- 
ample Q] Its identifying vector is v(X) = 1011000 and its 
extended representation is given by 

/ 1 1 1 \ 

pvt/ y\ _ [ 1 1 1 

1 J "0010101 

Voooioii/ 



The extended representation is redundant since the RREF 
define a unique subspace. Nevertheless, this representation will 
lead to more efficient enumerative coding. Some insight for this 
will be the following well known equality given in JT] p. 329]. 

Lemma 2: For all integers q, fc, and n, such that fc < n we 
have 

(4) 
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' n - 1 
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The order defined in Section [TTT] is based on Lemma [2] (applied 
recursively). Note that the number of subspaces in which 



1 is 



and the number of subspaces in which 



v{X) 1 = is q k 

Remark 3: A simple connection between OJ and (01 was 
given in (29] p. 68]. 

III. Coding Based on Extended Representation 

In this section we define a lexicographic order for the 
Grassmannian based on the extended representation. We present 
an enumerative coding technique for the Grassmannian using 
this order and discuss its complexity. 



Lemma 3: For 1 < j < n we have 



N 



Vj . 






n-j 


Xj . 




)= 


k - Wj 



Proof: Let X be a fc-dimensional subspace in G q (n, k) for 
which the first j columns in the extended representation are 
given by ^ ^ Xi ) ' Then in the last n — j entries of 

v(X) there are k — Wj ones, and the wj last rows of n — j last 
columns of EXT(X) have only zeroes. Therefore, restriction 
of EXT(X) to the first (k + 1) — Wj rows of the last n — j 
columns defines a subspace in Q q (n — j,k — Wj). Hence, we 
have 



N 



Xt 



n- j 
k — Wi 



A. Order for Q q {n, k) Based on the Extended Representation 

Let {x} denote the value of x — {x\,x<z, ...,x r ) £ 7L r q (or 
x = (x±,X2, x r ) T £ Zg), where the vector x is viewed as a 
number in base-q notation. Let {i} q be the base-q representation 
of the nonnegative integer i. The resulting vector is either a row 
vector or a column vector depending on the context. 

Let X,Y £ Gq(n, k) be two fc-dimensional subspaces and 
EXT(X), EXT(y) be the extended representations of X and 
Y, respectively. Let i be the least index such that EXT(X) 
and EXT(F) have different columns. We say that X < Y if 
{ } < { V } • Clearly, this definition induces an 

order for Q q {n, fc). 

Example 3: For X, Y, Z £ 2 (6,3) whose EXT(X), 
EXT(F) and EXT(Z) are given by 



EXT(X) 



EXT(y) 



EXT(Z) 



we have Y < X < Z. 



/111000X 
1 1 
1 
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B. Enumerative Coding Based on Extended Representation 

Let N ^ ^ be the number of elements in 

G q (n, fc) for which the first j columns in the extended rep- 
resentation are given by ^ ^ "' ) • 

Remark 4: We view all the q-ary vectors of length fc + 1 as 
our finite alphabet. Let S be the set of all q-ary (fc + 1) x n 
matrices which form extended representations of some fc- 
dimensional subspaces. Now, we can use Cover's method to 
encode/decode the Grassmannian. In this setting note that 
N ( " ^ is equivalent to ns{x\, X2, Xj), where 

^ has the role of Xi. 

Let Wj denotes the weight of the first j entries of v(X), i.e., 



Theorem 3: Let X £ G q {n, fc) be a subspace, where 
EXT(X) = 



x n 



v 2 «i 
Xi X\ 



Then the lexicographic index (decoding) of X, I EXT (JT), is given 
by 



i EXT{ x)^Y,(m k - W] - i +(i~v 3 )^f 1 1 ' 

3 = 1 



n-j 
k — Wj — i 



(5) 



Proof: By (0 we have that I EXT (^) is equal to 

N( I ^ - »' |. (6, 



3=1 1 « 



To compute the jth summand of (O, we distinguish between 
two cases. 

Case 1: Vj = 1. It implies that Xj has weight one, and 
its bottom Wj-i + 1 entries (as a column vector) are an 
one followed by uij-i zeroes, i.e., Xj = {q Wj ^ 1 } q - Hence, 
EXT(X) has the form 

V n ... Vj + 1 1 Wj-i ... Vx 

X n ... X, . , {</" : |, Xj.x ... Xi 

Therefore, a subspace Y £ G q (n, k) is lexicographically pre- 
ceding X, where EXT(y) has the same first j — 1 columns as 
EXT(X), if and only if EXT(F) has the form 

v' n ... v' j+1 Vj-i ... vi 
Y n ... Yj + i Yj Xj-i ... Xi 

Note, that Yj has zeroes in the last entries (since the 

leading coefficients of the last Wj-i rows are contained in 
(Xj_i • • • Xi)). The first fc — Wj-\ entries of Yj can have 
any values. 

Therefore, in this case the jth summand of © is equal to 



"3-1 -1 



E 



N 







Vj-l 



which is equal by Lemma [3] to 



Xi 



n-j 



k — w 



Case 2: Vj — 0. Since w 



3-1 



3-1 



— Yst=i ^ follows that the 



(7) 



4 



last Wj-i entries of Xj are zeroes, i.e., {Xj} is a multiple of Theorem 4: Encoding Algorithm A finds the subspace X E 



Hence, EXT(X) has the form 



t'n 

x n 



V 3 + l 

x, 



v 



3-1 



Xj Xj 



Xi 



Therefore, a subspace Y E G q (n, k) is lexicographically pre- 
ceding X, where EXT(y) has the same first j — 1 columns as 
EXTpf), if and only if EXT(Y) has the form 



Y„ 



where < s < 



U 3+l 







Y j+1 {s-q w ^} q Xj-x 



1. 



Thus, in this case the jth summand of (O is equal to 

q n( 



s=0 



s-q w ^} q X 



3-1 



which is equal by Lemma [3] to 



n - j 

k — Wj-i 



(8) 



Finally, combining equations ^ and ^ in Case 1 and Case 2 
implies equation ©. ■ 



Example 4: Let X G Q% (6, 3) be a subspace represented by 

EXT(X) 



/ 1 1 1 \ 

110 1 
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V 1 1 / 



By Theorem [3] we have that 



+ 2 2 



3 
2 

928. 



Now, suppose that an index < i < 



is given. En- 



coding Algorithm A finds X E G q (n, k) such that l EXT (X) = i. 

Encoding Algorithm A: 
Set io = i, wq = 0. 
For j = 1, 2, n do 

• if Wj-i = k then set vj = v(X)j = 0, Wj — Wj—i, 
Xj~ {0} q , and ij = ij-i; 

• otherwise 



- if i 



> 



n-3 
k — Wj — i 



then set 



= v(X)j = 1, Wj = Wj-x + 1, Xj = {q w ^} 



and in 



k—Wn 



-3 ~ - ? 

- otherwise let val = 

Vj = v(X)j = 0, w 
and ij = ij—i 
Form the output 



n-3 
k — Wj— 



ij-i/ 



k — w j _ i 



and set 



j - Wj-i, Xj = {val*q w i- 1 } Q , 
val * 



n-j 
k — w j _i 



EXT(X) 



W 2 Ul 
^2 Al 



Q q (n,k), such that I EXT (X) 

Proof: First we will show that the output of the algorithm 
is a fc-dimensional subspace. In other words, we will prove 
that the weight w n of identifying vector of the resulting 
subspace X is equal to k. We observe that the first "if of the 
algorithm implies that w n < k. Note also that ij > for all 
1 < j < n - Suppose that w n = k — t for some t > 0. Let 
n — k + 1 < j' < n be the last index where v(X)j> = 0. Then 
Wj' = k — t — n + j' = Wji—\. According to the algorithm, 

[ ] r ^'[r + ^3'i=° 

(since t > 0), which contradicts the observation that ij > 
for each 1 < j < n. 

Let 5j be the jth summand of I EXT (X), given in Q, i.e., 
l EXT (X) = Yjt=i ^ P rove me theorem it is sufficient to 
show that ij = i — J2t=i &t f° r all 1 < j < tt. and i n = 0. The 
proof will be inductive. 

By the algorithm, for each coordinate 1 < j < n — k, 



ij_ x - q 



k—Wj—i 



{Xj} 



n-3 
k — Wj—i 

n-j 
k — Wj—i 



, ifv(X) 3 =l 
, if v(X)j = 



Thus. 



lj_X 



v(X) jq 



k — Wj 



n-j 
k — Wj-i 



(1 - v(X)j) 



{X3} 



n-j 

k — Wj—l 



= ij.x ~ Sj (9) 



for all 1 < j < n — k. Thus, for j = 1 we have ix 



Sx. 



We assume that i n - 



Et=i S t> ^ 3 > !• By <EK 



■Sj+i, therefore, ij +1 — i— Yjt=i St~S^ 



■j+X 



3+1 - i ~J2lii &t- 

Now, we will show that for all < j < n, ij is the 
lexicographic index of a subspace in Q q {n — j,k — Wj) with 
given j first columns of its representation matrix. It will 
complete the proof since i n is the index of subspace in G q (0, 0) 
and thus it is equal to 0. 

for all < 



It is sufficient to prove that ij < 

j < n. The proof will be inductive. For j = we observe that 

' -3+1 ' 
- «y-i 



i < 



is given. Assume that ij-x < 



We will show that L < 



n-3 



two cases. 
Case 1: ij—x > 
algorithm, Vj 

ij-x - g fc - tu ^ 1 
n-j + l ] ,, 

k — w j _ 1 

< 



k - 
„k—Wj-i 
1, Wj 

n - j 



We distinguish between 



n-3 
k — wj _j 

= W 3-l 



Then, by the 

1, and ij = 



By the assumption, 



< 



n-3 
Wj-l 



- w 



3 

i-l 
3 

'j . 



Case 2:ij-x < q k ~ w ^- 1 



n-3 



k - ■ 



and thus by Lemma [2] 



, Then, by the algorithm, 



0, Wj 



Wj 



lj-X - 



x, and 



lj-XI 



k - ■ 



y-i 



k - 



-3 
"3-1 



5 



ij-il 


n-j 






n-j 




n-j 


k — w j _ i 






k — Wj — i 


q 


k — Wj—i 



since we can write |_§J < a < ([fj + 1)6 for all positive 
integers a and b. ■ 

Example 5: Let q = 2, n — 6, fc = 3, and i — 928. By 
using the Encoding Algorithm A we will find the subspace 
X € £2(6,3) such that I EXT (X) = i. We apply the following 
steps of the algorithm. 

j = l:i = 928 < 2 3 



7 val = L928/155J 
153. 

j = 2: h = 153 



5 

5,Xi 



1240 and hence vx = v(X)i 
, andii = 928-5-155 



> 



u 2 = v(X) 2 = 1, X 2 = 
j = 3: i 2 = 33 > 2 2 
X 3 = f 1 ) , and i 3 
j = 4: i 3 = 5 < 2 1 
m/ = [5/3J = 1, X 4 
j = 5: i 4 = 2 > 2 1 

^« = (s 

j = 6: w 5 



= 120 and hence 

and i 2 = 153 - 120 = 33. 
28 and hence V3 = v(X)s 
33 - 28 = 5. 

1 = 6 and hence vx = v(X)x 
1 J 2 

= ( J, and i 4 = 5 - 3 = 2. 

[ = 2 and hence i; 5 = v(X) 5 
L J 2 

and « 5 = 2 - 2 = 0. 
3 = fc and hence i>6 = v(X)q = 0, 
X 6 = ( ^ , and i 6 = i 5 = 0. 

Therefore, we obtain a subspace X E £2(6,3) whose 
extended representation is given by 



1, 



= 0, 



= 1, 



EXT(X) 



/ 1 1 1 

110 1 

1 

V 1 1 



C. Complexity 

We consider the complexity of computation of lexicographic 
index I EXT (-) in (0. Note that all the integers that we use in the 
calculations are q-ary integers. Let M[a, b] denotes the number 
of operations for the multiplication of two q-ary integers of 
length a and b. It is known QUI p. 634], that for a > b, 
M[a, b] = a log b log log b. 

First, we calculate the length of the g-ary integer which 
represents the largest Gaussian coefficient in (O. This Gaussian 
coefficient is 

_ (q 71 - 1 - 1) • • • (q n - k - 1) 



n- 1 

k 



(q k -l)---(q- 1) 



and hence this length is less than k(n — k). 
If Wj = Wj—i then 



n-j 




' n - (j + 1) " 




k — Wj-i 


1 


k — Wj 





i—k—j-\-Wj 



If Wj = Wj- 



1 then 



n-j 




' n - (j + 1) " 




k — Wj — \ 




k — Wj 






1 




Q 



-.it—I, 



- 1 



(10) 



(ii) 



The Gaussian coefficients in (f5J) can be derived from the 
identifying vector. Their computation is done by ( fTOb and ( fTTT >. 
Hence, the complexity for computation of all the Gaussian 
coefficients that we need in © is 0(nM[k(n — k),n}). 



Since multiplication or division by q 1 is done by a shift of i 
digits, there are n — k indices where Vj = 0, and the length of 
{Xj} is k, it follows that the complexity of these operations is 
0((n — k)M[k(n — k), k]). Finally, in (0 there are at most n 
additions of integers whose length is at most fe(n — k+ 1), and 
therefore the complexity of these operations can be omitted. 

Hence, the complexity of computation of I EXT (-) in (0 is 
0(nM[k(n — k),n]), i.e., 0(nk(n — k) log n log log n). 

Therefore, we have proved the following theorem: 

Theorem 5: The computation complexity of the lexico- 
graphic index (decoding) in (|5]l is 0(nk(n — k) log n log log n) 
digits operations. 

If k < log n log log n then the Gaussian coefficients in (0 
can be computed more efficiently. For their computation we can 
use Lemma [2] To compute , we need to compute 71 

L -I flf l K 1 q 

for all -q and k such that < k < k and 0<r/ — k < n — k. 
It requires at most k(n — k) additions of integers whose length 
is at most k(n — k), and a total of at most k(n — k) shifts. All 
other computations do not change and can be omitted from the 
total complexity. Thus, we have 

Theorem 6: If mm{k,n — k} < log n log log n, then the 
computation complexity of the lexicographic index in (0 is 
0(n 2 min{fc, n — fc} 2 ) digits operations. 

Finally, in a similar way we can show that the computation 
complexity of Encoding Algorithm A is the same as the 
computation complexity given for the decoding in Theorem 
and in Theorem [6] 

IV. Coding Based on Ferrers Tableaux Form 

In this section we present an enumerative coding for the 
Grassmannian based on the Ferrers tableaux form represen- 
tation of fc-dimensional subspaces. Note, that even so this 
enumerative coding is less efficient, it is more intuitive and 
might have its own applications. Lexicodes based on the related 
order, were found to be larger than the known codes 



A. Enumerative Coding for Ferrers Diagrams of the Same Size 

Let be a Ferrers diagram of size m embedded in a k x (ri- 
fe) box. We represent T by an integer vector of length n — fc, 
(J- n -k, J~2, where Ti is equal to the number of dots in 
the i-th column of J 7 , 1 < i < n — fc. Note, that the columns 
are numbered from right to left and that < J-^+i < T% < fc 
for alll<i<n — k — 1. Let T and T be two Ferrers diagrams 
of the same size. We say that T < T if Ti > Ti for the least 
index i such that T\ 7^ Ti, i.e., in the least column where they 
have a different number of dots, T has more dots than T . This 
is similar to the lexicographic order defined in the literature for 
unrestricted partitions, e.g. iFJTI . Il32l pp. 93-98]. 

Let 7V m (J r J , J-~2, Fx) be the number of Ferrers diagrams 
of size m embedded in a fc x (n — fc) box, for which the first 
j columns are given by (Tj, F^Fi). 

Lemma 4: If 1 < j < n — fc and < m < k(n — fc) then 

3 

Nm^j,..., ^2,^1) =p(Fj,n- fc- j,m- J^Fi). 

i=l 

Proof: The lemma is an immediate consequence from 
the fact that T = (F n -k, J- 2 , Fx) is a Ferrers diagram 
with m dots embedded in a fc x (n — fc) box if and only if 
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(T n -k, —, Fj+i) is a Ferrers diagram with m — J2l=i T dots 
embedded in an Tj x (n — k — j) box. ■ 

Remark 5: We view the set Z^+i = {0,1,..., k} as our 
finite alphabet since < T < k. Let S be the set of all (n—k)- 
tuples over Zk+i which represent Ferrers diagrams embedded 
in a k x (n — k) box. In other words, (J- n -k, ■■■,J'2,J'i) £ S 
if and only if < T < i < k for each 2 < i < n — k. 
Now, we can use Cover's method to encode/decode the set of 
Ferrers diagrams with to dots embedded in a k X (n — k) box. 
In this setting note that N m (Tj , . . ., T2 , T\) is equivalent to 
ns(xi, X2, ■ ■ ■ , Xj), where T has the role of x*. 

Theorem 7: Let T = (T n -k, J~2, T) be a Ferrers dia- 
gram of size to embedded in a k x (n — k) box. Then the 
lexicographic index (decoding), ind m , of T among all the 
Ferrers diagrams with the same size m is given by 

ri — fe Fj — i j — 1 

ind m (J") = ^2 E P( a : n ~ ^-j^-^T- a), (12) 

j — 1 a— J-j+1 i=l 

where we define To = k. 

Proof: By (O we have that 

n— k Fj—i 

ind m (J") = ^2 E Nm ( a ' ^i-ii • 77 2' -^i)- 
3=i 0=^+1 

The theorem follows now from Lemma |4] ■ 
Remark 6: The summation in Theorem [7] is over larger 
values, while the summation in (0 is over smaller values, due 
to the defined order (T < J- if Ti > Ti for the least index i). 

Theorem |7] implies that if we can calculate p(k,r),m) effi- 
ciently then we can calculate ind^T 7 ) efficiently for a Ferrers 
diagram of size m embedded in a k X (n — k) box. 

Now suppose that an index < i < p(k, n — fc, to) is given. 
Encoding Algorithm B finds a Ferrers diagram T of size to 
embedded in a k x (n — k) box, such that ind m (J r ) = i. 

Encoding Algorithm B: 

Step 1: Set J-'q — k, l\ = 0, h = i, iq = i; 

. while h > N m (T„ - t x ) set h = h - N m (T - ti), 
h=h + i; 

• set T\ = J-o — l\, and i\ = h; 

Step 2: For j = 2, n — k do 

• ^ = m men set = 0; 

• otherwise do 
begin 

- set £j = 0, h = 

- while h > Nm^Tj-i — £j, J^j—i, set h = h — 
N m (Tj-i - i 3 ,Tj-i, J 7 !), tj = ij + l; 

- set Tj = J-j-i — ij, and ij = h; 
end {begin} 

Step 3: Form the output T = (F n -k: J-2, 

Remark 7: We did not join Step 1 and Step 2, since 
N m {Fj-i — lj,Tj-\, ■■■,J r i) is not defined for j = 1, 

B. Order for Q q {n, k) Based on Ferrers Tableaux Form 

Let X, Y 6 fc) be two fc-dimensional subspaces 

and let Tx, be the related Ferrers diagrams. Let x = 

(x 1 ,x 2 ,-.,x\jr xl ) and y = (y lt j/ 2 , V\r Y |) be the entries 



vectors of .F(X) and T(Y), respectively. These entries are 
numbered from right to left, and from top to bottom. 

We say that X < Y if one of the following conditions holds. 

. \J r x\> |7V|; 

. \T X \ = \T Y \ and T x < Ty\ 
• Tx = Ty and {x} < {y}. 
Clearly, this definition induces an order for Q q (n, k). 
Example 6: Let X, Y, Z,W £ £2(6, 3) be given by 



F(X) 



HZ) 



111 
111 
1 



111 
1 1 




T(Y) 



10 1 



1 1 



111 
1 1 
1 



Tz = T w and by definition Z < W. Clearly, \Tx\ = \Ty\ > 
\T Z \ and Ty < T x . Thus, Y < X < Z < W. 

C. Enumerative Coding Based on Ferrers Tableaux Form 

In this subsection, we use the given order of Ferrers tableaux 
forms and Theorem [2] for enumerative coding for Q q (n,k). 

Theorem 8: Let X £ G q (n, k), Tx be the Ferrers diagram 
of X, and let x = (xi, x 2 , %\f x \) be the entries vector 
of T(X). Then the lexicographic index (decoding) of X, 
Indjr(X), defined by the order based on Ferrers tableaux form, 
is given by 



Ind^(X) = 



k(n— k) 

E 

■-\F X \+1 



+md { r x \(T x )q^ +{x}, (13) 



where a^, \T X \ + 1 < i < k(n — k), is defined in Theorem [2] 
Proof: To find Indjr(X) we have to calculate the number 
of fc-dimensional subspaces which are preceding X according 
to the order defined above. 

1) All the fc-dimensional subspaces with Ferrers diagrams 
which have more dots than Tx are preceding X. Their 
number is 2Ji=\F x \+i «»? • 

2) There are ind|jr x | (Tx) Ferrers diagrams with \Tx\ 
dots which are preceding X. Hence, there are 
indijr x | {T x )q' :Fx ' fc-dimensional subspaces whose Fer- 
rers diagrams have \Tx\ dots and preceding X. 

3) Finally, the number of fc-dimensional subspaces whose 
Ferrers diagram is Tx which are preceding X is {x}. 

■ 

Example 7: Let X £ ^2(6, 3) be the subspace of Example[U 
whose Ferrers tableaux form and Ferrers diagram are 



1 1 

and Tx 
1 



T(X) = 

By Theorem [8] we have that 

9 

lndr(X) - a ^ + ind 4 (^)2 4 + {(1011)}. 

i=5 

Since a 5 = 3, a 6 = 3, 07 = 2, as = 1, ag = 1 (see |Q] pp. 
326-328]), md 4 (T x ) = 0, and {(1011)} = 11, it follows that 
Indjr(X) = 1323. 

Now suppose that an index < i < ? is given. 

Encoding Algorithm C finds a subspace X £ Q q \n, k) such 
that Indjr(X) = i. 
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Encoding Algorithm C: 
Set io = i. 

For j = 0, . . . , k(n — fc) do 
• if u < 



k(n-k)-j 



then 



set 



\Fx\= k(n- k) - j, F x 



ind 



( L „k(n-k)-j J )> 



{zj— [ fc( „!ffc)_j J g fc ^" " is assigned to x (the entries 

vector of F(X)) and stop; 

otherwise set = i i - a k ( n -k)-j<l ki - n ~ k ^ : ' ■ 



D. Complexity 

We consider the complexity of the calculation of the lex- 
icographic index lndjr(X), for X £ Q q (n,k), whose Ferrers 
diagram is Fx — (Fn-k, We will use the following 

lemma concerning partitions to find a bound on the length of 
q-ary integers which represent the value of p(k, n — k, i). 

Lemma 5: For any given n, k, and i, we have 
p(k, n — fc, i) < e^V^. 

Proof: Clearly, p(k,n — k,i) < p(i), where p(i) is the 
number of unrestricted partitions of i. It is known [1, p. 160] 
that p(i) < e^v'T* and the lemma follows. ■ 

Theorem 9: The computation complexity of the lexico- 
graphic index (decoding) in ( TBI is 0(fc 5 / 2 (n — fc) 5 / 2 ) digit 
operations. 

Proof: First, we combine the expressions in (fT2l i and (fT~3t 
to obtain: 



Ind^(X) 



k(n—k) 

^2 p(k,n 



k,i)q l + {x} 



-k T 3 



J2 J2 P{a,n-k- j,\F x \-^^~a). (U) 



j=l a=J r j + l 



i=l 



By the recurrence relation of Lemma Q] we can compute the 
table of p(j, i, i) for j < k, £ < ry, and i < m with no 
more than mkrj additions. By Lemma [5] each integer in such 
addition has 0(y/k(n — fc)) digits. Therefore, the computation 
of all the values which are needed from the table takes 
0(fc 5 / 2 (n — fc) 5 / 2 ) digit operations. 

The number of additions in (TBI is 0(k(n— k)). Each integer 
in this addition has 0(k(n — fc)) digits (as a consequence of 
Lemma [5] and the powers of q in (TT4b). The multiplication by 
q l is a shift by i symbols. Hence, these additions and shifts do 
not increase the complexity. ■ 

Similarly, we can prove the following theorem. 

Theorem 10: The computation complexity of Encoding Al- 
gorithm C is 0(fc 5 / 2 (n — fc) 5 / 2 ) digit operations. 

Remark 8: If k(n — fc) — \Fx\ is a small integer then the 
complexity of the computation becomes much smaller than 
the complexity given in Theorems [9] and [10] For example, 
if \Fx\ = k(n — fc) then the complexity of the enumerative 
decoding is 0(k(n — fc)) since Indjr(X) = {x} in (TBI . 

It is worth to mention in this context that the number of 
operations in the algorithms can be made smaller if we will 
consider the following two observations ll27l p. 47]: 

• If mi < m2 < then p(k, rj, mi) < p(k, 77, m-i). 

• p(k,rj ) rn) = p(k,r),kr] — m) and hence we can assume 
that m < !f. 



V. Combination of the Coding Techniques 

By Theorems |5J [6] and |£l it is clear that the enumerative 
coding based on the extended representation is more efficient 
than the one based on Ferrers tableaux form. But, for some of 
fc— dimensional subspaces of F™ the enumerative coding based 
on Ferrers tableaux form is more efficient than the one based 
on the extended representation (see Remark [HJ. This is the 
motivation for combining the two methods. 

The only disadvantage of the Ferrers tableaux form coding 
is the computation of the c^'s and indi^.1 (Fx) in Theorem[8] 
This is the reason for its relatively higher complexity. The 
advantage of this coding is that once the values of the a^'s 
and the value of indi^i (Fx) are known, the computation of 
Indjr(X), for X g Q q (n, fc), is immediate. Our solutions for 
the computation of the Qj's and \ad.\j^ x \(Fx) are relatively not 
efficient and this is the main reason why we suggested to use 
the enumerative coding based of the RREF and the identifying 
vector of a subspace. The only disadvantage of this enumerative 
coding is the computation of the Gaussian coefficients in (0. It 
appears that a combination of the two methods is more efficient 
than the efficiency of each one separately. The complexity will 
remain 0(nk(n — fc) log n log log n), but the constant will be 
considerably reduced on the average. This can be done if there 
won't be any need for the computation of the a^'s and the 
computation of ind^i^A') will be efficient. 



It was proved in [9] that q 



k(n—k) 



< 



< 4q k ( n ~ k *> for 



< fc < n. Thus, more than | of the fc-dimensional subspaces 
in G q (n, fc) have the unique Ferrers diagram with k(n — fc) 
dots, where the identifying vector consists of fc ones followed 
by n — fc zeroes. All the codewords of the Reed-Solomon-like 
code in [9] have this Ferrers diagram. Note that most of the 
fc-dimensional subspaces have Ferrers diagrams with a large 
number of dots. We will encode/decode these subspaces by 
the Ferrers tableaux form coding and the other subspaces by 
the extended representation coding. We will choose a set Sjf 
with a small number of Ferrers diagrams. Sjr will contain the 
largest Ferrers diagrams. The Ferrers tableaux form coding will 
be applied on these diagrams. 

We say that a subspace X g Gq(n, fc) is of Type Sjr if Fx € 
Sjf- In the new order these subspaces are ordered first, and their 
internal order is defined as the order of the Ferrers tableaux 
forms in Section [IV] The order of the other subspaces is defined 
by the order of the extended representation in Section [HI] We 
define a new index function I com b as follows: 



comb 



Ind^(X) 



l EXT (X) + Ax(Sjr) otherwise ' 



(15) 



where Ax(«SV) is the number of subspaces of Type which 
are lexicographically succeeding X by the extended represen- 
tation ordering. These Ax(Sj^) subspaces are preceding X in 
the ordering induced by combining the two coding methods. 

We demonstrate the method for the simple case where Sjr 
consists of the unique Ferrers diagram with k(n — fc) dots. 

Lemma 6: Let Sjr be a set of Ferrers diagrams, embedded 
in a fc x (n — fc) box, which contains only one Ferrers diagram, 
the unique one with k(n~ fc) dots. Let X G G q (n, fc), X ^ S-f, 
RE(X) = (X n , . . . ,Xi), and let I, < £ < n - k - 1, be 
the number of consecutive zeroes before the first one (from 
the right) in the identifying vector v(X). Then Ax (Sjr-) = 

ElivV-i-Wk 



k(n — k — i) 



s 



Proof: If I — then v(X)i = 1 and hence there are no 
subspaces of Type Sjr which are lexicographically succeeding 
X and hence A x (Sj) = 0. For 1 < I < n - k - 1, let 
X±, Xg be the first £ columns of RE(X). All the subspaces 
of Type Sj- in which the value of the first column is greater 
than {^fi}, are lexicographically succeeding X. There are 
(q k — 1 — {Xi})q k ( n ~ k ~ 1 ^ such subspaces. All the subspaces 
of Type Sjr in which the first i — 1 columns, 2 < i < n — k — 1, 
are equal to the first i — 1 columns of RE(X), and the value 
of the zth column is greater than {X}, are lexicographically 
succeeding X. There are (q k — 1 — {X;})g fe ("~ fc-J ) such 
subspaces. Therefore, there are Yli=i(l k ~ 1 — {Xi})q k ^ n ~ k ~'^ 
subspaces of Type Sj which are lexicographically succeeding 
X by the extended representation ordering. ■ 
Example 8: Let X be the subspace of Example [4] By 
Example [4] we have Iex T (X) = 928, and by Lemma [6] we 
have Ax(Sjr) = (2 3 - 1 - 5)2 3 2 = 2 7 . Hence, I com bP0 = 
WPO + A x (5^) = 928 + 128 = 1056. 



is given. 



Now, suppose that an index < i < 

Based on (fl3T l and Lemma [6] we can find the subspace X such 
that Icomb(^) = i, where Sj- consists of the unique Ferrers 
diagram with k(n — k) dots. We omit the details of the encoding 
algorithm. 

VI. Conclusion 

Three methods of enumerative coding for the Grassmannian 
are presented. The first is based on the representation of 
subspaces by their identifying vector and their reduced row 
echelon form. The second is based on the Ferrers tableaux form 
representation of subspaces. The complexity of the first method 
is superior on the complexity of the second one. The third 
method is a combination of the first two. On average it reduces 
the constant in the first term of the complexity compared to the 
complexity of the first method. Improving on these methods is 
a problem for future research. 

The enumerative coding is based on an order for the Grass- 
mannian related to a specific representation. This order can be 
used to form lexicographic codes Il25ll in the Grassmannian. 
To our surprise some of these lexicographic codes form the 
best known error-correcting codes in the Grassmannian. For 
example, a lexicode of size 4605 in £2(8,4) with minimum 
subspace distance 4 (see [9| for the distance definition) was 
generated based on Ferrers tableaux form order (compared to 
the largest previously known code of size 4573 generated by a 
multilevel construction [18]). These codes also revealed a new 
method to form error-correcting codes in the Grassmannian. 
This topic is considered in ll26l . 

Construction of a lexicode might require to generate all 
subspaces of Q q (n, k) by the given lexicographic order. Usually, 
this does not require to use the enumerative coding since the 
subspaces are generated one after another. By using one of 
our orders it is not difficult to prove that given a subspace 
X G G q [n, k), it takes no more than 0(kn) digit operations to 
generate the next subspace. 

Acknowledgment 

We thank the anonymous reviewers whose comments have 
helped to improve the presentation of this paper. 



[10 
[11 

[12 

[13 

[14 

[15 

[16 

[17 
[18 

[19 

[20 
[21 
[22 

[23 

[24 
[25 

[26 

[27 

[28 
[29 

[30 

[31 

[32 



References 

J. H. van Lint and R. M. Wilson, A course in Combinatorics, Cambridge 
University Press, 2001 (second edition). 

D. E. Knuth, "Subspaces, subsets, and partitions ," /. Combin. Theory, 
vol.10, pp. 178-180, 1971. 

S. Thomas, "Designs over finite fields", Geometriae Dedicata, vol. 21, 
pp. 237-242, 1987. 

W. J. Martin and X. J. Zhu, "Anticodes for the Grassman and biliniar 
forms graphs", Designs, Codes, and Cryptography, vol. 6, pp. 73-79, 
1995. 

S. Thomas, "Designs and partial geometries over finite fields," Geome- 
triae Dedicata, vol. 63, pp. 247-253, 1996. 

R. Ahlswede, H. K. Aydinian, and L. H. Khachatrian, "On perfect codes 
and related concepts," Designs, Codes, Crypt., vol. 22, 221-237, 2001. 
M. Schwartz, T. Etzion, "Codes and anticodes in the Grassman graph," 
Journal of Combinatorial Theory, Series A, vol. 97, pp. 27^12, 2002. 
M. Braun, A. Kerber, and R. Laue, "Systematic construction of g-analogs 
oft — (v, k, A)-designs," Designs, Codes, and Cryptography, vol. 34, pp. 
55-70, 2005. 

R. Koetter and F. R. Kschischang, "Coding for errors and erasures in ran- 
dom network coding," IEEE Trans. Inform. Theory, vol. 54, no. 8, 
pp. 3579-3591, August 2008. 

S. T. Xia and F. W. Fu, "Johnson type bounds on constant dimension 
codes," Designs, Codes, and Cryptography, vol.50, pp. 163-172, 2009. 
T. Etzion and A. Vardy, "Error-correcting codes in projective space", proc. 
Int. Symp. on Inform. Theory, Toronto, pp. 871-875, July 2008. 

F. Manganiello, E. Gorla, and J. Rosenthal, "Spread codes and spread 
decoding in network coding", proc. of Int. Symp. on Inform. Theory, pp. 
881-885, July 2008. 

D. Silva, F. R. Kschischang, and R. Koetter, "A rank-metric approach to 
error control in random network coding," IEEE Trans, on Inform. Theory, 
vol. IT-54, pp. 3951-3967, September 2008. 

D. Silva and F. R. Kschischang, "On metric for error correction in 
network coding," IEEE Trans. Inform. Theory, vol. IT-55, pp. 5479-5490, 
December 2009. 

M. Gadouleau and Z. Yan, "Constant-rank codes and their connection to 
constant-dimension codes," IEEE Trans. Inform. Theory, vol. IT-56, pp. 
3207-3216, July 2010. 

M. Gadouleau and Z. Yan, "On the decoder error probability of bounded 

rank distance decoders for maximum rank distance codes," IEEE Trans. 

on Inform. Theory, vol. IT-54, pp. 3202-3206, July 2008. 

M. Gadouleau and Z. Yan, "Construction and covering properties of 

constant-dimension codes," arxiv.org/abs/0903.2675 . 

T. Etzion and N. Silberstein, "Error-correcting codes in projective space 

via rank-metric codes and Ferrers diagrams", IEEE Trans. Inform. Theory, 

vol. IT-55, pp. 2909-2919, July 2009. 

A. Kohnert and S. Kurz, "Construction of large constant dimension codes 
with a prescribed minimum distance," Lecture Notes Computer Science, 
Vol. 5393, pp. 31^12, 2008. 

V. Skachek, "Recursive code construction for random networks," IEEE 
Trans. Inform. Theory, vol. IT-56, pp. 1378-1382, March 2010. 
T. M. Cover, "Enumerative source encoding," IEEE Trans. Inform. The- 
ory, vol.IT-19, no. 1, pp. 73-77, Jan. 1973. 

V. Braun and K. A. S. Immink, "An Enumerative coding technique for 
DC-free runlength-limited sequences," IEEE Trans, on Commun., vol. IT- 
48, no. 1, pp. 2024-2031, Dec. 2000. 

O. f. Kurmaev, "Enumerative coding for constant-weight binary sequences 
with constrained run-length of zeros," Problems of Inform. Tran., vol. 38, 
no. 1, pp. 249-254, 2002. 

K. A. S. Immink, Codes for Mass Data Storage Systems, Shannon 
Foundation Publishers, 1999. 

J. H. Conway and N. J. A. Sloane, "Lexicographic codes: error-correcting 
codes from game theory," IEEE Trans. Inform. Theory, vol. IT-32, pp. 
337-348, May 1986. 

N. Silberstein and T. Etzion, "Large constant dimension codes and 
lexicodes," arxiv.org/abs/1003.4879.2010. 

G. E. Andrews, The Theory of Partitions, Cambridge University Press, 
1984. 

R. P. Stanley, Enumerative Combinatorics, volume I, Wadsworth, 1986. 
G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge University 
Press, 2004. 

D. E. Knuth, The Art of Computer Programming, Vol.2, Seminumerical 
Algorithms, Third Ed., Addison- Wesley, 1997. 

T. V. Narayana, R. M. Mathsen, and J. Sarangi, "An algorithm for 
generating partitions and its applications," J. Combinatorial Theory, vol. 
11, pp. 54-61, 1971. 

F. Ruskey, Combinatorial Generation, Working Version, University of 
Victoria, Victoria, Canada, 2001. 



Natalia Silberstein was born in Novosibirsk, Russia, in 
1977. She received the B.A. and M.Sc. degrees from the 
Technion - Israel Institute of Technology, Haifa, Israel, in 2004 
and 2007, respectively, from the Computer Science Department 
and the Applied Mathematics Department, respectively. She is 
currently working toward the Ph.D. degree in the department 
of Computer Science at the Technion. Her research interests 
include algebraic error-correction coding, coding theory, and 
combinatorial designs. 



Tuvi Etzion (M'89-SM'94-F'04) was born in Tel Aviv, 
Israel, in 1956. He received the B.A., M.Sc, and D.Sc. degrees 
from the Technion - Israel Institute of Technology, Haifa, Israel, 
in 1980, 1982, and 1984, respectively. 

From 1984 he held a position in the department of Computer 
Science at the Technion, where he has a Professor position. 
During the years 1986-1987 he was Visiting Research Professor 
with the Department of Electrical Engineering - Systems at 
the University of Southern California, Los Angeles. During 
the summers of 1990 and 1991 he was visiting Bellcore in 
Morristown, New Jersey. During the years 1994-1996 he was a 
Visiting Research Fellow in the Computer Science Department 
at Royal Holloway College, Egham, England. He also had 
several visits to the Coordinated Science Laboratory at Univer- 
sity of Illinois in Urbana-Champaign during the years 1995- 
1998, two visits to HP Bristol during the summers of 1996, 
2000, a few visits to the department of Electrical Engineering, 
University of California at San Diego during the years 2000- 
2010, and several visits to the Mathematics department at Royal 
Holloway College, Egham, England, during the years 2007- 
2009. 

His research interests include applications of discrete mathe- 
matics to problems in computer science and information theory, 
coding theory, and combinatorial designs. 

Dr Etzion was an Associate Editor for Coding Theory for 
the IEEE Transactions on Information Theory from 2006 till 
2009. 



